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Abstract 

The atomic and magnetic structures of (CuX)LaNb207 (X=C1 and Br) are investigated 
using the density- functional calculations. Among several dozens of examined structures, an 
orthorhombic distorted 2x2 structure, in which the displacement pattern of X halogens 
resembles the model conjectured previously based on the empirical information is identified as 
the most stable one. The displacements of X halogens, together with those of Cu ions, result 
in the formation of X-C\x-X-Cn-X zigzag chains in the two materials. Detailed analyses of 
the atomic structures predict that (CuX)LaNb207 crystallizes in the space group Pbam. The 
nearest-neighbor interactions within the zigzag chains are determined to be antiferromagnetic 
(AFM) for (CuCl)LaNb 2 7 but ferromagnetic (FM) for (CuBr)LaNb 2 7 . On the other hand, the 
first two neighboring interactions between the Cu cations from adjacent chains are found to be 
AFM and FM respectively for both compounds. The magnitudes of all these in-plane exchange 
couplings in (CuBr)LaNb207 are evaluated to be about three times those in (CuCl)LaNb207. In 
addition, a sizable AFM inter-plane interaction is found between the Cu ions separated by two 
NbOg octahedra. The fourth-neighbor interactions are also discussed. The present study strongly 
suggests the necessity to go beyond the square J\ — J2 model in order to correctly account for the 
magnetic property of (CuX)LaNb207. 
PACS: 71.15.Mb, 75.45.+j 
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I. INTRODUCTION 



Low- dimensional quantum spin systems with frustrated interactions have drawn consid- 
erable attention for several decades 1|. In particular, the square-lattice 5 = 1/2 frustrated 
Heisenberg magnets with first-neighbor exchange constant J\ and second-neighbor constant 
Jo are increasingly interesting due to their unusual ground states and quantum phenomena 
[2]- 13}. Based on the J\ — J 2 model studies, there exist several phases as a function of 
J2/J1. When J 1 dominates or J 2 is ferromagnetic (FM), the system is either Neel anti- 



13|). when 



ferromagnetic (NAFM) or FM depending on the sign of J\ (Refs. 
J 2 is antiferromagnetic (AFM) and dominates, there appears the so-called columnar AFM 
(CAFM) order [3|, |7J with antiferromagnetically coupled FM chains. The CAFM and FM, or 
CAFM and NAFM ordered phases are separated by the intermediate quantum-disordered 
phases, the nature of which is not yet fully resolved {2], U47, 12, 13|. The recent discoveries 



of quasi-two-dimensional materials are realizations to test the validity of the J\ — J 2 model. 
Prominent among them are Li 2 VO(Si,Ge)04 (Refs. 
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16|-|l8j), (CuBrlA'^a 



Zr, Sr, Ba) (Refs. 
and (CuX)LaNb 2 7 (X=C1, Br) (Refs. 
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15J), AB(VO)(P0 4 ) 2 (A,S=Pb, 
A'=Ca, Sr, Ba, Pb; £'=Nb, Ta) (Ref. (lj), 
2l|). (CuX)LaNb 2 7 compounds are of par- 



ticular interest because they allow systematic tuning and understanding of the structural 
and magnetic properties, which are plausibly connected with the phenomenon of high-T c 
superconducting cuprates. 

Although divalent copper with the electronic configuration d 9 should be Jahn- Teller active 
and lead to the cooperative lattice distortion (e.g., perovskite KCuF 3 (Ref. (22 1 ) ) , the precise 
crystal structure of the layered copper oxyhalides (CuX)LaNb 2 7 is still under debate. 
Earlier structural studies on (CuX)LaNb 2 7 were carried out with the tetragonal space 



group P^/mmm, where the Cu and X sites possess the C4 symmetry [23|, |24j. While the 
Rietveld refinement gave satisfactory results, the thermal parameter for halogens remained 
large. Besides, in this structure copper is in a significantly squeezed octahedral coordination 
with two short Cu-0 bonds (about 1.9 A) and four rather long Cu-X bonds (2.7 A), which 
are also quite unusual. Subsequently, the neutron diffraction experiment [25j proposed 
that the CI ions in (CuCl)LaNb 2 7 ((CuCl)LNO) shifted away from the ideal Wyckoff 16 
position [23]. The transmission electron microscopy measurement on (CuCl)LNO (Ref. 2p| ) 
revealed superlattice reflections corresponding to an enlarged 2x2 unit cell. The nuclear 
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magnetic resonance and the nuclear quadrupole resonance experiments for (CuCl)LNO and 



(CuBr)LaNb 2 7 ((CuBr)LNO) 
at both Cu and CI/Br sites [20 



urther demonstrated the lack of the tetragonal symmetry 



2l|. 



The magnetic properties of (CuCl)LNO and (CuBr)LNO are also unusual and lack a 



clear microscopic interpretation. 
, Q that are incompatible 
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he former exhibits a spin liquid phase with a spin gap 



it has been reported 



28j with the square J\ — J2 model. On the other hand, 



29] that the replacement of CI by Br leads to a CAFM order in 



(CuBr)LNO at low temperatures. However, it is unclear whether the Cu ions connected 
with the dominant exchange interaction couple ferromagnetically or antiferromagnetically 



21]. Moreover, both (CuCl)LNO and (CuBr)LNO are claimed to be FM J x compounds 



29] whose justifications largely rely on the J\ — J 2 model. Yet, the structural study [21] 



raised doubts over the validity of the model. Therefore, unambiguous determination of the 
crystal structure is crucial for understanding these complex systems. 

At present, there are several structural models proposed for the CuX plane. Whangbo 
and Dai [28] suggested a model that consists of different ring clusters to explore the ex- 
change couplings. However, the existence of inequivalent Cu and CI sites in such a model 
is in contradiction to the experimental results that both Cu and X occupy a unique crys- 
tallographic site with no substantial disorder 20j, |21]. Yoshida et al. [20], based on the 
empirical evidence, proposed an orthorhombic distorted 2x2 structure (hereafter referred 
to as the YY model). In this model, the displacement of CI ions generates different exchange 
couplings among the nearest neighboring Cu pairs. A Cu dimer formed by the dominant ex- 



change interaction was considered 
model was shown 



to study the spin-gap behavior. The same structural 



21] to consistently account for (CuBr)LNO. The third model, suggested 
by Tsirlin and Rosner (TR) [3^] is also characterized by an ordering pattern but with a 
2x1 periodicity, where the local environment of copper is distorted to form the CUO2CI2 
plaquette. 

First-principles calculations have proven to be an appealing method to deal with complex 
systems 3ll433]. Such a method can efficiently and reliably calculate the total energy, which 
is crucial in determining the most stable structure in order to study all relevant physical 
properties. In this work, we will investigate the atomic structure and resultant magnetic 
property of (CuCl)LNO and (CuBr)LNO based on the density functional theory. Our results 
show that, among several dozens of examined structures, the distortion pattern of the most 
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stable one is similar to that of the YY model. The displacement of the A ions changes 
the environment of copper to form the CUO2A2 plaquette. In addition, these two materials 
crystallize in the space group Pbam. The FM chains in CAFM (CuBrjLNO are found to 
be along the direction which is contrary to the previous conjecture [21] . It will be shown 
that (CuCl)LNO still belongs to the AFM J\ compound. The first- and second-neighbor 
exchange couplings of (CuCl)LNO and (CuBr)LNO are also discussed in detail. 



II. CRYSTAL STRUCTURE AND COMPUTATIONAL DETAILS 



Figure [T] illustrates the basic crystal structure of the copper oxyhalides (CuA)LaNb20y. 
It is made up of copper-halogen planes and nonmagnetic double-perovskite LaNb207 slabs. 
The La ions are located at the 12-coordinate sites of the double-perovskite slabs. The CuA 
planes and the LaNh^Oy slabs alternate along the c direction such that the copper is six- fold 
coordinated, bridging between the apical O ions of NbC>6 octahedra and surrounded by four 
X halogens. Because of the short Cu-0 bond length (~1.9 A), the CuA plane is more 
appropriately considered as a CUAO2 layer. The initial structural study on (CuA)LaNb20y 
was carried out with the space group P^/mmm, where both Cu and X have the C4 symmetry 
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24j . In this model (hereafter referred to as C4), the Cu and X ions are located at the 



Wyckoff Id and 16 positions, respectively (Fig. 0a)). Later studies [20|, |25] proposed that 
CI ions are displaced from the ^-symmetry positions. The YY 2x2 model is represented 
in Fig. [2](b). The displacement of X ions on the CuA plane leads to the formation of 
the A-Cu-A-Cu-A zigzag chains, as indicated in Fig. [2](b). The original equivalent and 
perpendicular Cu chains are now distinguishable. Here, the direction extending along the 
zigzag chains is defined as the b axis. 

The present calculations were based on the generalized gradient approximation (GGA) 



34| to 



;he exchange-correlation energy functional of the density functional theory. It is 



known 



35M37II that Cu-derived oxide compounds are usually strongly correlated systems. 



The correlation effect is important for the present systems to understand their ground state. 



Therefore, the on-site Coulomb interaction U for Cu 3d electrons was also included 
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(GGA+tf) in this work. Since the on-site exchange interaction J is expected to be less 



influenced by the solid state effects 



30j . the relation J = 0.1U was used 



for different 



choices of U. The projector-augmented-wave potentials, as implemented in VASP 
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4Q|, 
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were employed for the interactions between the ions and valence electrons. The plane-wave 
basis set with an energy cut-off of 500 eV was used. To minimize numerical uncertainties, 
structural optimizations were performed using a 2 x 2 supercell for all the test structures 
unless specified otherwise. The 6x6x4 Monkhorst-Pack grids were taken to sample the 
corresponding Brillouin zone. The lattice parameters and atomic positions were relaxed 
until the total energy changed by less than 10 -6 eV per conventional cell and the residual 
force was smaller than 0.01 eV/A. 



III. RESULTS AND DISCUSSION 



A. Energetics 



We first calculated the total energies of the C4 structure and several 2x1 and 2x2 
distorted structures with different displacement patterns of the X halogens. The YY 2x2 
model is found to be the most stable one. As compared to the C4 structure, the YY model 
has a significant 0.3 and 0.2 eV/fu lowering in the energy of (CuCl)LNO and (CuBr)LNO, 
respectively. This directly rules out the possibility that the two compounds crystallize in 
the C4 symmetry. Particularly, over the full U range from eV to 8 eV, the YY model is 
0.1 eV/fu lower than the TR 2x1 model (see Fig. 4 in Ref. 30] ) for both materials. The 
nonmagnetic calculations with and without U lead to the same conclusion. Therefore, the 
structural distortion outweights the magnetism and on-site correlation effects in determining 
the atomic structure. Note that, besides Refs. and 21], very recent experimental 

evidences 41] also confirm that the original unit cell should be double along both the a and 



b axes for the family of these compounds. Our study therefore provided theoretical support 
for the stabilization of the YY 2x2 model. 

To examine whether there exists other more stable structure with the X ions restricted 
to the Cu plane, we perform the calculations for (CuCl)LNO with twenty sets of random 
displacements of all four CI ions from the positions in the YY model. However, no such 
structure was found. The resultant configurations of the trial structures are either relaxed 
back to the YY model or trapped into a nearby higher energy minimum. 

Next, we allow the halogens in the YY model to move off the plane. It is found that the 
X ions in the relaxed structures are 0.02 - 0.04 A away from the Cu plane. However, the 
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change in the total energy is rather small. At U =0 and 8 eV, the results of both materials 
show that the energy differences are only within 2 meV per 2x2 supercell while the energies 
for the structures with the X ions fixed in the plane remain lower. We also examine the 
two structures in Figs. 16(e) and 16(f) of Ref. 20], which are based on another 2x2 
configuration with the CI ions displaced away from the Cu plane in a different way. The 
calculations indicate that, after relaxation, both are energetically about 0.2 eV/fu higher 
than the YY model. The increase in the total energy mostly comes from the different in- 
plane structural distortions. Again, the contribution from the ^-component shift of CI ions 
is rather minor. Hence, the distortion on the CuX plane is predominantly crucial to stabilize 



the atomic structure. In t 
ions kept in the Cu plane 



re following discussion, we shall focus on the YY model with X 



42| . 



Now, we analyze the total energies influenced by the on-site Coulomb interaction and 
the different magnetic configurations shown in Fig. [31 Here, SCI and SC3 are FM and 
NAFM. SC2 and SC4 are both CAFM, with the FM chain along the b and a directions, 
respectively. The results are displayed in Fig. HI where the energy of SC2 was chosen 
as a reference. For (CuCl)LNO, the energies of SCI, SC2, and SC4 are very competing. 
The differences among them are within 1 meV/fu when U > 6. The first two are even 
almost identical around U=4 eV. Clearly, Fig. Hta) shows that the SC3 is the lowest energy 
state and its energy is well separated from those of the other three magnetic structures. 
In the (CuBr)LNO case, similar tiny energy differences but between SCI, SC2, and SC3 
are also found. Interestingly, when the FM chain in the CAFM state is set parallel along 
the a axis, as in SC4, the total energy over the examined U range is much higher than 
those of the other three configurations, indicating that the Cu ions along the a-axis should 



not couple ferromagnetically. This finding is contrary to the previous conjecture [21(. The 
different energy ordering for the four magnetic configurations between the two compounds 
are conceivable since the magnetic interactions through the path Cu-X-Cvl depend subtly 
on the small structural variation via the X-ion size effect. We will return to this issue in 
Sec. HID when considering the various exchange couplings. 
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B. Atomic structure 



Tables [Jlists the fully optimized structural parameters of both materials. For comparison, 
those obtained by the C4 model are also included. As can be seen in this table, the evaluated 
lattice constants are in good agreement with the experimental data 24|, |25|. The discrepan- 
cies between them are only within 1%, the typical errors in the density-functional calculation. 
The a and b lattice constants of (CuCl)LNO are smaller than those of (CuBr)LNO, which 
is due primarily to the size effect of Br in the layered structure. 

To discuss the structural distortion, we take the (CuCl)LNO case as an example. In 
the C4 model, copper is in the squeezed octahedral coordination with four long Cu-Cl 
bonds [d(Cu-Cl)=2.77 A] and two short Cu-0 bonds [d(Cu-0)=1.85 A]. The displacements 
of the CI and Cu ions in the YY model lead to two shorter Cu-Cl bonds of 2.38 and 
2.39 A, forming the CI- Cu-Cl- Cu-Cl zigzag chain to stabilize the structure. The rest two 
Cu-Cl interatomic distances are increased to 3.27 — 3.29 A. In particular, the Cu-0 bond 
length remains short after the structural distortion (from 1.85 A to 1.88 A), indicating the 
strong bonding character between Cu and O ions. The calculated interatomic distances are 
comparable to those reported previously [251 ] . As a result, the distortion yields the nearly 
planar CUO2CI2 rather than the octahedral CUO2CI4 environment around the Cu ion (Fig. 
Ufa)). The resultant CUO2CI2 planar structure is reminiscent of the conventional CUO4 
plaquette, which is commonly observed in copper oxides, e.g., La 2 Cu0 4 (Ref. [36|) and 
Sr 2 Cu03 (Ref. [37|). It should be noted that the Cu02Cl2-plaquette zigzag chains was also 
reported in the TR model 30] . Additionally, the basic electronic structure is similar to that 
of the CUO4 planar unit, which will be demonstrated in the next section. Combined with 
the energetic advantage mentioned above, the YY model provides a realistic description for 
the atomic structures of (CuCl)LNO and (CuBr)LNO . 

From a closer analysis of the positions of all ions in (CuCl)LNO, we found that the 
distorted atomic structure in the YY model belongs to the space group Pbam (No. 55) 44]. 
The atomic positions are summarized in Table [III Clearly, the deviations of the CI ions from 
the C^-symmetry positions are as large as 0.66 A, and these values are four times larger 
than those of Cu ions. Note that the displacements along the a axis are more significant 
than those along the b axis for both ions to form the zigzag chains. 

As expected, the structural distortion on the CuCl plane leads La, Nb and O ions to 
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shift from the the ^-symmetry positions. Figure [5]^b) and Table [TT] show the significant 
tilting and distortion of NbOg octahedra. Such a tilting distortion is typical for perovskite 
oxides structures 45]. Particularly, La ions shift along the b axis by an amount of 0.10 A. 



This displacement of La from the ^-symmetry positions well agrees with the experimen- 
tal nonzero value of the EFG tensor at La sites [20], a strong evidence for the structural 
distortion in (CuCl)LNO. We also found that Nb ions shift along the a axis by a relatively 
smaller amount of 0.02 A. Note that the a (b) component of La (Nb) displacement is almost 
negligible. 

Taking into account the tilting of the NbOg octahedra in the (CuCl)LNO is important 
for providing a realistic description of the distortion on the CuCl plane. Figure El^c), the 
top view of the atomic structure, clearly demonstrates that the cooperative tilting of the 
NbC>6 octahedra in the space group Pbam results in a 2 x 2 periodicity and leads to the 
zigzag chains with the same periodicity. It is worth pointing out that the higher symmetric 
2x1 zigzag chains in the TR model were investigated without consideration of the effect 
due to the tilting of the Nb0 6 octahedra, where these octahedra were still kept at the 
C4 tetragonal sites [30]. Allowing the tilting distortion of NbOg octahedra in (CuCl)LNO 
lowers the symmetry of the atomic structure and correspondingly that of the zigzag chains 
and therefore leads to a lower total energy. In the YY model, the zigzag chains have the 
glide symmetry about u = 1/4 and v = 1/4. Specifically, the Cu-Cl bond of 2.39 A in the 
zigzag chain is next to the Cl-Cu bond of 2.38 A in the adjacent chain and vice versa. As 
compared to those in the TR model, such a complementary arrangement between adjacent 
chains in the YY model allows a further lowering in energy. Now, it is evident that the YY 
2x2 model in the present study is energetically more stable than the TR 2x1 one. Similar 
conclusion holds for (CuBr)LNO. 



C. electronic structure 



Figure [6] depicts the orbital- and site-projected densities of states (DOS) of (CuCl)LNO 
with U=6 eV, where the valence-band maximum (E v ) is set to zero. The orbitals are 
projected in the local coordinates with the x and y axes directed to the neighboring CI ions 
and the z axis coinciding with the crystal c axis (Fig. Elb)). Among the major valence 
state region of 6.3 eV, the higher-energy part consists almost exclusively of O and CI p 
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states. There is larger contribution from the CI p state just below E v . The lower-energy 
part, dominated by the Cu d states, is splitted into the doubly occupied dxy, dyz, dzx, 
d(x 2 — y 2 ), and singly occupied d(3z 2 — r 2 ) states. As compared to the GGA DOS (not 
shown here), the GGA+[7 shows an essential redistribution of the Cu 3d DOS, i.e., from 
being above to below the O and CI p states. That the energy gap (E g ) lies between the 
occupied anion p states and the unoccupied Cu d states is similar to those in the charge- 
transfer insulators, e.g., La2Cu04 (Ref. |36|) and Sr 2 Cu03 (Ref. |37|). The sharp peak of 
the low-lying Cu d(3z 2 — r 2 ) state is a result of the strong bonding between the Cu and O 
ions with a considerably short Cu-0 bond length of 1.88 A (see Tabled]). Note that the 
d(3z 2 — r 2 ) orbital was hybridized with little CI p component. The on-site correlation U leads 
to the half-filling of the d(3z 2 — r 2 ) orbital, of which the lobes point to the O ions. These 
results imply a single orbital ground state. Figure [6] shows that, due to the hybridization 
with the O p state j^f]], the d(3z 2 — r 2 ) bonding-antibonding separation (8.3 eV) is larger 
than the value of U. Therefore, the electronic structure due to the CUO2CI2 plaquette in 



the YY model is very similar to those of other copper oxides [36J with planar CUO4 units. 

We found that the structural distortion and magnetism together already open up the gap. 
The Eg of (CuCl)LNO obtained by the GGA is 0.3 eV. However, this result is insufficient 
to account for the green color appearance 24| of this compound. At U = 6 eV, the E g is 
increased to 1.8 eV. Further increase of U makes no significant widening for the band gap. 
The main structures in the DOS of (CuCl)LNO are also found in that of (CuBr)LNO, except 
for the smaller E g of 1.5 eV. At this choice of U, the local magnetic moment at the Cu site 
CuBr)LNO is evaluated to be 0.6 Bohr magneton, which agrees with the experiments 
29J. An amount of 0.1 Bohr magneton at Br sites is also observed. Hence, we choose 



of 
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the optimal U = 6 eV case to discuss the corresponding atomic and electronic properties. 



D. exchange interaction 

Finally, we discuss the exchange couplings for both (CuCl)LNO and (CuBr)LNO. In 
the undistorted C4 structure, the interactions between the Cu ions can be approximately 
modeled by the Heisenberg Hamiltonian H — J\ J2nn &i ' $j + Ji J22NN & ' ^jj where the 
sums run over the first and second nearest-neighbor pairs, respectively, and Si corresponds 
to the spin moment at site i. The relevant exchange couplings can be then determined by 
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applying the model to the calculated energies of different spin configurations. For the YY 
model, the formation of the X-Cu-X-Cu-X zigzag chains along the b axis (Fig. EJ^b)) lifts 
the tetragonal symmetry and leads to inequivalent superexchange pathways, as indicated 
in Fig. [2J As a result, the J\ in the C4 structure is split into Jn, J12, and J 13 , with the 
former two now being the first neighboring inter-chain interactions and the latter the first 
neighboring intra-chain interaction. The original J2 coupling is split into two inequivalent 
J 2\ and J22, which are correspondingly the second neighboring inter-chain interactions. We 
investigate these interactions via the various spin configurations in Fig. |3j The results are 
summarized in Table IHIl 

We first discuss the results from the C4 model. Table HTT1 shows that J\ is almost negligible 
as compared to J 2 - This is reasonable because, as illustrated in Fig. [7(a), there is no overlap 
between the Cu2 d(x 2 — y 2 ) and X4 p orbitals. Therefore, even with the obvious overlapping 
of the Cul d and X4 p orbitals, Cul and Cu2 could hardly interact with each other. On the 
other hand, Cul can interact with Cu3 via the X4 p orbital. 

Based on the J\ — J2 model, both (CuCl)LNO and (CuBr)LNO were previously claimed 
271 . l29j ] to be FM Ji magnets with competing AFM J2 interactions, as in the case of 
Pb 2 VO(P0 4 ) 2 (Ref. [18|). Table DID indeed shows that J x < and J 2 > for both materials 
in the C4 model, a direct consequence of the Hund's coupling and virtual electron hopping. 
However, the recent structural study [21] has raised serious doubt over the validity of the 
J\ — J2 model in such materials. Our calculations also indicate that consideration of the 
structural distortion leads to the opposite results. For (CuCl)LNO, the effective interactions 
( Jn + J12 + 2 Ji 3 ) /4 and ( J 2 i + J22) /2 in the YY model are found to be AFM and FM, respec- 
tively. And they both become FM for (CuBr)LNO at large Us. These results come from the 
complicated interplay between the Hund's coupling, virtual electron hopping, the distorted 



structure and X-ion size effect. It should be noted that the TR model [30J] also results 
in a leading AFM coupling in (CuCl)LNO. This implies that the simple J\ — J2 model is 
unable to describe the present systems. Moreover, the first neighboring interactions become 
more significant as compared to the second neighboring ones. Figure W[b) clearly shows 
that, unlike the C4 case, L Cul-A4-Cu2 is no longer 90° (Table HJ due to the structural 
distortion. This will lead to the overlapping of Cu2 (Cu4) d and X4 p orbitals, and enhance 
the interaction between Cul and Cu2 (Cu4). 

In fact, Fig. |^a) shows that for (CuCl)LNO, (J u + J 12 )/2 > 0, J 13 > 0, and 
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(J21 + ^22)/2 < for all the Us considered. It is now clear that, since the interactions 
due to all the corresponding spin pairs in SC3 satisfy these conditions, (CuCl)LNO in SC3 
is much more stable than in the rest three configurations of Fig. [3J Actually, SC3 is the 
most stable structure among all the spin configurations with the interactions up to second- 
nearest neighbors. However, by comparing Fig. |S(b) with Fig. E(a), we found that the 
J13 in (CuBr)LNO becomes FM. This could be ascribed to the fact that the interaction J13 
depends sensitively on the angle of the Cul-A4-Cu2 superexchange path in Fig. and the 
replacement of CI by Br will change this angle and modify the interaction. Therefore, for 
(CuBr)LNO, the first neighboring couplings within the chains are FM, and the first two 
neighboring couplings between adjacent chains are AFM and FM, respectively. None of the 
four structures in Fig. [3] satisfies these conditions. Specifically, the interactions due to the 
corresponding spin pairs in SC4 are all opposite to these couplings, giving rise to the result 
of the SC4 being the highest-energy structure for (CuBr)LNO. 



The above analysis seems to indicate that, contrary to previous expectations [20 
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27|, 



29[], (CuCl)LNO rather than (CuBr)LNO is less frustrated [29j. To check the reliability of our 



calculations, we perform the total-energy calculation for an additional structure SC5 with 
three of the four spins being the same but opposite to the fourth one. For all the possible 
choices of four magnetic structures containing SC5 in solving the Heisenberg Hamiltonian, 
the deviations of the relevant couplings (dashed lines in Fig. [8]) from those obtained by 
SCI— SC4 are less than 1.0 meV. More importantly, the signs of these interactions remain 
unaltered for both materials. 

One plausible explanation for the above puzzling discrepancy is that the third nei ghb oring 
couplings between different zigzag chains 46j may not be completely negligible 27], 130 ] . 
Actually, in the YY model, the two CUO2A2 plaquettes with the Cul and Cu5 ions in Fig. 
W[h) could be considered approximately coplanar. Kageyama et al. 47j argue that this kind 
of coplanarity provides an opportunity for the interaction between Cul and Cu5 through 
the overlap of the Cul d(x 2 — y 2 ) — X4 p — X5 p— Cu5 d(x 2 — y 2 ) orbitals. So, from the 
structural geometry point of view, such a long-range coupling (8.5 A) could be possible. 
However, to examine and identify these couplings, one has to take into account additional 
eight inequivalent couplings and use a larger supercell whose corresponding calculations are 
very time-consuming and yet, likely, not accurate enough for the present purposes (Fig. Sj). 
Therefore, we will not address this issue presently. Further theoretical and experimental 
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work is required to clarify this point. 

It is worth pointing out that the first neighboring interactions within the chains and 
between the adjacent chains have the opposite signs for (CuBr)LNO. The spatial asymmetry 
of these results again demonstrates the inappropriateness of the square J\ — J2 model for 
the Br compound. Furthermore, for both materials, the couplings between the adjacent 
chains are very competitive to the intra-chain interactions, sharply contrary to the previous 



conclusion 



301 ] . Tsirlin and Rosner have argued 30] that in the TR model, where the basic 



structure element is also the CuC^CVplaquette zigzag chain, the large hopping runs along 
the chain and leads to the strongest interaction. According to their discussion, the inter- 
chain interaction is rather weak due to the long "nonbonding" Cu-Cl distance and the lack 
of the proper superexchange path. In the present study, despite the similar backbone in 
the YY model, the couplings between the adjacent chains are shown to be still substantial. 
As mentioned before, the strength of exchange interactions between two spin sites should 
be determined by the overlap of orbitals rather than the distance between them. The 
interactions between the Cu ions from adjacent chains could be significant through the path 
mediated by the extended 3p orbital of CI ions (vs 2 ~) and would be enhanced in the Br 
case with the further extended 4p orbital. Indeed, our calculations show that all the in-plane 
exchange couplings in (CuBr)LNO are three times larger than those in (CuCl)LNO. 

For the inter-layer interaction Jj_, Table HTT1 shows that the J± is AFM, in agreement with 
the experiment [291. When compared to the in-plane interaction, the J± in (CuCl)LNO is 
non- negligible, implying that some long-path (12 A) interaction between the Cu ions is still 
cooperative. The origin of this long-range coupling could be associated with the interaction 
through the Cu d(3z 2 — r 2 ) orbital. As discussed in Sec. IIIC, this orbital strongly overlaps 
with the O p z orbital. The O p z orbitals further couple with Nb 4d orbitals. Therefore, 
the inter-plane coupling Jj_ shall involve the Cu-O-Nb-O-Nb-O-Cu path. In (CuBr)LNO, 
however, the coupling is found to be relatively less significant. The strength of all the 
interactions interested here is decreased with increasing U. The evolution is expected since 
adding U makes the wavefunctions more localized and the virtual electrons hopping less 
favorable. 
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IV. CONCLUSIONS 



In conclusion, we have investigated the atomic structure and magnetic property of the 
copper oxyhalides (CuY)LaNb 2 07(X=Cl and Br) based on the density functional theory. 
The calculations show that, among the examined structures, the YY 2x2 model proposed 



by Yoshida et al. |20j has the lowest energy. This model is significantly more stable than 
both the undistorted C4 structure and the TR 2x1 model suggested recently by Tsirlin 
and Rosner 30(. The X and Cu ions in the YY model are displaced to form the X-Cu- 
X-Cu-X zigzag chains and the local environment of copper is distorted to form a nearly 
Cu0 2 X 2 plaquette. We found that (CuX)LaNb 2 7 crystallizes in the space group Pbam. 
The cooperative tilting of the NbOg octahedra leads to the lower symmetry of the zigzag 
chains with a 2 x 2 periodicity. With consideration of the on-site Coulomb interaction, the 
YY model shows the single-orbital scenario typical for copper oxides and oxyhalides. 

We concluded that (CuCl)LNO is still the AFM J\ magnet with mixing FM J 2 inter- 
actions. For (CuCl)LNO, the first neighboring interactions within the zigzag chains are 
AFM, and the first two neighboring couplings between adjacent chains are AFM and FM, 
respectively. However, the replacement of CI by Br modifies the first neighboring intra-chain 
interaction to be FM for (CuBr)LNO. Despite the "well" -separated zigzag chains in the YY 
model, the couplings between adjacent chains are comparable to those within the chain. The 
opposite signs of the inter- and intra-chain interactions in (CuBr)LNO reflect the spatial 
asymmetry and therefore the failure of the simple J\ — J 2 model for such material. All the 
in-plane exchange couplings in (CuBr)LNO are shown to be three times those in the CI 
counterpart. It is found that the inter-plane interaction J± is AFM, in agreement with the 
experiment [29]. The present study strongly suggests that the simple square J\ — J 2 model 
should be modified to explore the magnetic property of CuXLaNb 2 07. We hope the present 
calculations will shed light on the precise crystallographic determination and the magnetic 
properties of CuXLaNb 2 07. 
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FIGURE CAPTIONS 



Fig. 1: (Color online) Crystal structure of (CuX)LaNb207 (X=C1, Br) in the tetragonal 
space group P^/mmm. 



Fig. 2: (Color online) (a) The undistorted tetragonal model and (b) the distorted 
model proposed by Yoshida et al. 2CJ for (CuX)LaNb207 (X=C1, Br). Large and small 
spheres denote X and Cu ions, respectively. The relevant exchange couplings are also 
indicated. 



Fig. 3: (Color online) The four different spin configurations of Cu ions considered 
in the present study. Large and small spheres denote CI/Br and Cu ions, respectively. 



Fig. 4: (Color online) Total energies of (a) (CuCl)LaNb207 and (b) (CuBr)LaNb20y in the 
spin configurations shown in Fig. [31 The energies are relative to that of SC2. 



Fig. 5: (Color online) Perspective view of (a) the Cu02X 2 -plaquette zigzag chains and (b) 
tilted Nb0 6 octahedra of (CuX)LaNb 2 7 (X=C1, Br) in the space group Pbam. (c) Top 
view of (b). Here, the symbols for the various atomic species are the same as those in Fig. HJ 



Fig. 6: (Color online) Orbital- and site-projected density of states (DOS) of (CuCl)LaNb207, 
obtained by U — 6 eV and SC2 in Fig. [3J The energy is relative to the valence-band 
maximum. 



Fig. 7: (Color online) Schematic plot of Cu d(x 2 — y 2 ) and X p orbitals of (CuX)LaNb 2 07 
(X=C1, Br) in the space group (a) P^/mmm and (b) Pbam. 



Fig. 8: (Color online) The ^/-dependence of the in-plane exchange couplings of (a) 
(CuCl)LaNb 2 7 and (b) (CuBr)LaNb 2 7 . The solid lines are calculated from SC1-SC4 in 
Fig. [3j The dashed lines show the uncertainty of the calculations. See the text for details. 
Note that the scale of exchange coupling in (a) is only half of that in (b). 
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TABLE I: The calculated lattice constants (A), relevant interatomic distances (A), and bond angles 
(°) of (CuX)LaNb 2 7 (X=Cl, Br), obtained by U = 6 eV and SC2 in Fig. El The latter two refer 
to those in Fig. [2(b). In the bottom part, the results in the undistorted tetragonal model are listed 
for comparison. The values in parentheses are the corresponding experimental data. 





(CuCl)LaNb 2 7 


(CuBr)LaNb 2 7 


a 


7.868 


7.889 


b 


7.883 


7.914 


c 


11.878 


11.853 


Cul-Xl 


2.39 (2.40t) 


2.54 


Cul-X2 


3.27 (3.14+) 


3.11 


Cul-X3 


3.29 


3.11 


Cul-X4 


2.38 


2.52 


Cul-0 


1.88 (1.84+) 


1.88 


Z Xl-Cul-X2 


83.6 


81.5 


Z X2-Cul-X4 


87.1 


88.8 


Z X4-Cul-X3 


102.7 


101.3 


Z X3-Cul-Xl 


86.6 


88.3 


Z Cul-X4-Cu2 


112.2 


103.4 


a* 


3.914 (3.884+) 


3.942 (3.899++) 


c* 


11.892 (11.736+) 


11.853 (11.706++) 


Cu-X* 


2.77 


2.79 


Cu-O* 


1.85 


1.87 



n 

f: Ref. [25 
tf: Ref. M 
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TABLE II: The calculated atomic structural parameters of (CuX)LaNb207 (X=C1 and Br) in 
the space group Pbam, obtained by U = 6 eV and SC2 in Fig. u, v, and w denote fractional 
coordinates based on the a, b, and c lattice constants, respectively. 



(CuCl)LaNb 2 7 (CuBr)LaNb 2 7 



ion 


site 


u 


V 


w 


u 


V 


w 


Cu 


Ah 


0.2706 


0.0077 


0.5 


0.2720 


0.0057 


0.5 


X 


Ah 


0.4185 


0.2711 


0.5 


0.4481 


0.2713 


0.5 


La 


4g 


0.5000 


0.2622 





0.5000 


0.2624 





Nb 


8i 


0.2522 


0.0000 


0.1911 


0.2520 


0.9998 


0.1903 


01 


4/ 


0.5 





0.1336 


0.5 





0.1334 


02 


4e 








0.1841 








0.1834 


03 


8? 


0.2498 


0.2501 


0.1522 


0.2498 


0.2499 


0.1520 


04 


4ff 


0.2004 


0.0001 





0.2004 


0.0000 





05 


8? 


0.2830 


0.0008 


0.3417 


0.2818 


0.0000 


0.3413 
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TABLE III: The exchange couplings (meV) of (CuCl)LaNb207 and (CuBr)LaNb207. The notation 
is explained in the text. U (eV) is the on-site Coulomb correlation interaction. 







(CuCl)LaNb 2 7 








(CuBr)LaNb 2 7 








U (Jii 


+ ^12)/2 Jl3 {J21 


+ «/ 22 )/2 J± 


Ji J2 


(J11 + Ji 2 )/2 




(J21 + ^22)/2 J± 


Ji 


h 





8.9 


21.3 


-3.7 


5.1 


-1.8 18.4 


28.2 


-16.8 


-10.6 


5.4 


-2.9 


-1.5 


4 


3.6 


5.4 


-1.8 


2.0 


-3.6 26.9 


11.9 


-11.9 


-5.9 


2.2 


-1.8 


23.3 


6 


2.2 


3.0 


-1.0 


1.3 


-2.5 19.4 


7.9 


-8.6 


-3.5 


1.5 


-0.1 


25.5 


8 


1.7 


2.1 


-0.7 


0.7 


-1.7 13.5 


5.3 


-6.2 


-2.0 


0.9 


-0.3 


18.6 
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